The Z Transformation 

By H. A. HELM 

(Manuscripl received April 15, 1058) 

The Slieltjcs integral is used to develop a rigorous derivation of the z 
transform. Sufficient properties of the transformation are included to form 
a reasonably complete basis for the operational solution of constaid coeffi- 
cient, linear, finite difference equations. 

I. INTRODUCTION 

The desire to use digital computers in automatic control loops created 
the need for methods with which to analyze systems that are partly 
continuous and partly discrete. Since the methods of network theory 
could be applied to the analysis of the continuous part of such a hybrid 
system, it was natural that such methods should be extended to in- 
clude the discrete case. This resulted in the z transform introduced by 
Raggazini and Zadeh. 1 There is today an extensive literature devoted 
to the z transform."' 3 ' 4 However, the fundamental assumption of the z 
transform derivation is that the process of instantaneous sampling is 
equivalent to the amplitude modulation of a train of unit impulses by 
the "sampled" function. But the unit impulse as commonly defined has 
infinite height and zero width, and the process of amplitude modulating 
such a function is not intuitively clear. While it is true that such a 
process may be considered as an approximation to the behavior of a 
linear network with an amplifier and sampling switch, "impulse sam- 
pling" bears no simple relation to the manner in which the digital com- 
puter operates. The digital computer, in the type of real time operation 
typical of control system applications, works with sequences of numbers 
which represent a continuous function evaluated at particular instance 
of time. Since these numbers must of necessity be finite, "impulse sam- 
pling" is not an obvious mathematical model for describing the working 
of the computer. It is the intention of this paper to define the problem 
from the point of view of operations within the computer and to develop 
a rigorous and appealing derivation of the z transform. 

In place of impulse modulation, the alternate approach is taken 
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of generalizing the definition of the Laplace transformation by means of 
the Stieltjes integral. This approach has the advantage of rigor and of 
more closely relating the operational solutions of continuous and dis- 
crete systems. As developed below, only rational transforms are con- 
sidered. Also, in general, functions involving derivatives of impulses can- 
not be represented by the Stieltjes integral. Practically, these restrictions 
do not limit the applications greatly. Derivation of the principal proper- 
ties of the z transform, based upon the Stieltjes integral definition, are 
given in the Appendix. 

II. THE LAPLACE-STIBLTJES TRANSFORMATION 

The Stieltjes integral has the important property of including both 
sums and limits of sums (integrals). For the reader's convenience the 
definition of the Stieltjes integral as given in Widder 5 is repeated. 

Let an interval (a, b) be divided into sub-intervals in the following 
maimer: 

a = x < x\ < Xo < • • • < x n = b 

and let A equal the largest of these subintervals. Then the Stieltjes in- 
tegral of f(x) with respect to a(x) from a to b is 

f(x) dct{x) = lira Z/(n)[a(.r,) - aU-i)], (1) 



/ 



where Xk-i ^ f* ^ .r* . The left-hand side of (1) is the usual notation for 
a Stieltjes integral. The integral itself is defined only when the limit on 
the right exists. It can be shown" that the integral exists if f(x) is con- 
tinuous and if a(x) is monotonic but not necessarily continuous, i.e., 
nonin creasing or nondecreasing. We shall assume that both these condi- 
tions apply to all functions to be considered. However, these conditions 
are quite strong and will be somewhat relaxed subsequently. 

It is now possible to generalize the Laplace transformation by making 
the defining integral a Stieltjes integral. Thus, 

L.[f(t)] = I" f(t)e-' da(t), (2) 

Jo 

where 

S = a + ju. 

It is assumed that (2) is subject to all the above restrictions. As defined 
by (2) the Laplace-Stieltjes transformation actually defines a different 
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transformation for each different selection of the function a(t). lia(t) = t, 
(2) reduces to the usual definition of the Laplace transform. If a(l) is 
continuous and has a continuous derivative, (2) reduces to the Laplace 
transform of «'(/)/(/), which may he handled by the usual theorems of 
the Laplace transform. However, if a(t) is not continuous a new class 
of transformations results. For the purpose of this paper, the func- 
tion a(t) will be defined as the "staircase" function which increases by 
unity at integral multiples of T but remains constant between such 
points. This function is shown in Fig. 1. The constant T is equivalent 
to the sampling period of modulation theory. From this point on, the 
function a(t) will be assumed to be that of Fig. 1. Thus, (2) may be 
evaluated for this a(t) by means of (1) as 



f f{t)e-'Mt) = Z f(nT)e- nTl = /<V'), 

Jtt n =0 



(3) 



where a(t) is given by Fig. 1, s = a + j<a and f(nT) is the function /(/) 
evaluated at t = tiT. Since «(/) is monotonic, and we have restricted 
/(/) to continuous functions, (2) clearly exists. This does not imply that 
the series on the right converges or diverges. Since no part of /(/) be- 
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Fifl. 1 — "Staircase" function a (I). 
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tween multiples of T affects (3), f(t) need only be continuous and well 
defined in the neighborhood of the points t = nT; i.e., have no discon- 
tinuities or "jumps" at multiples of T. 

A simple change of variable in (3) introduces the z transform. Let 
z = e sT and substitute in (3) : 

Ls[f(t)] = f f(t)z-" T MO = E f(nT)z-' = f{z). (4) 

JO n =0 

The expression (4) emphasizes the power series nature of the z transform. 
Since the functions with which we shall deal are convergent, they can 
almost always be written in closed form. In fact, to allow the order of 
certain limit functions to be interchanged in the development of various 
theorems, absolute convergence of (4) will be assumed. That is, 

E \f(nT)z-\ =S M, (5) 

n=0 

where M is finite although possibly very large. The magnitude of z is 
( r " Ta , and if (5) holds for some a, say o- , it will obviously hold for any 
a > (To . In the following sections, the functions /(/) are now restricted to 
those functions for which (5) holds. Rewriting (5) as a Stieltjes integral 
leads to: 

[ | f(t) | (T al da(t) ^ M, (5a) 

Jo 

which holds for the Laplace transform when a(t) = /. The restriction of 
absolute convergence has little practical effect upon the utility of the 
transform. 

III. THE INVERSE TRANSFORM 

An inverse transform is necessary for operational completeness. The 
derivation is straightforward, making use of Cauchy's method for evalu- 
ating the coefficients of a power series. 7 The proof proceeds from the 
definition of the z transform in (4) : 

W(t)] = E f(nT) z -' = F{z). 

;i=0 

Expanding, 

f(z) = /(0) + .f(7>-' + • • • + f(nT - T)z- ( "- l) + f{nT)z~ n 

+ /(, i 7'+ 7>- < " +1 > + ••• , 
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which is absolutely convergent for | z | ^ e'° T . Since F(z) is a power 
series, absolute convergence also implies uniform convergence within 
the radius of convergence. Hence, (0) may he integrated term by term 
along a contour such that \z\ > e"° T . Multiplying by z"' 1 and so inte- 
grating gives: 

J z'- l F(z) dz = (/(0)* n_1 dz + [ r f(T)z n - 2 dz+ ■■■ 

+ ff{nT - T) dz + f /(nT)z- 1 dz + • ■ ■ (7) 

+ f f(nT + mT)z- (m+l) ■■■, 
Jv 

where m = 1, 2, 3, • • • and, Y the contour of integration, is a circle en- 
closing the origin of the z plane and whose radius is greater than e'° . 

It is obvious that all integrals except / J{nT)z' 1 dz either have no singu- 
larity within the contour of integration, and hence are zero, or else have 
a pole at the origin of order greater than unity, and hence also are zero. 
Therefore, (7) reduces to 

f z"- l b\z) dz = f(nT) f z' 1 dz, (8) 

since/(//7') is the function /(0 evaluated at I = nT, and hence a constant- 
Thus, the inversion formula becomes 



KnT) - -'. I z n ~ 1 F{z)dz. (9) 

'lirl h 



I 

It is necessary to add a word of caution here. Equation ( ( J) only gives 
the value of the function at one point, / = nT. Obviously, by assigning 
any given integer to n the value at any point may be obtained and, 
usually, this is sufficient. However, in the case of certain summations 
this will lead to a very confusing notation. To avoid this, we introduce 
the notation 

\f(nT)\ m 
to indicate the sequence 

/(()),/( T), f(2T) ■■■ J(nT), ■■■ J(mT). 

Omission of the subscript m will denote an infinite sequence. 
It is now possible to summarize the Laplace-Stieltjes pairs: 
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L,[f(t)]= [" J(t)er" da(t) = F(e T ) (3) 



or 

L.IM) = f f(t)z~" T Mt) = F{z), (3a) 

Jo 

W\F(z)] = ±-. f z- l F(z) dz = f{nT). (9) 

2tj Jr 

A closer inspection of the above pairs indicates a very interesting prop- 
erty of the Laplace-Stieltjes transformation. A continuous function f(l) 
is transformed and then the inverse operation performed. However, the 
function is then only defined at integral multiples of T, i.e., at / = nT. 
This is exactly equivalent to "sampling"; that is, the computer, by 
means of some encoding device, evaluates instantly a continuous func- 
tion of time (commonly represented by voltages, shaft positions, etc.) 
at periodic intervals. The L„ transformation is a mathematical model or 
abstraction which represents this process. It is very often the case that 
the function which is to be sampled (a voltage, for instance) is applied 
to a linear network and then sampled, and it would be very convenient 
to be able to analyze the system directly from the Laplace transform. 
Such results follow from the relationships between the Laplace and 
Laplace-Stieltjes transforms which we develop below. 

IV. THE RELATIONSHIP BETWIOUN THE LAPLACE AND LAPLACE-STIELTJES 
TRANSFORMATIONS 

The transform pairs (3) or (3a) and (9) are sufficient to derive opera- 
tional methods for discrete linear systems which are similar to those of 
continuous linear systems. However, there are several reasons for ex- 
amining the relationship between the two transforms. In the first place, 
many of the most interesting problems arise from the analysis of systems 
which are partly continuous and partly discrete. Also, the relationship 
between the Laplace-Stieltjes and Laplace transforms may be expressed 
as a convolution integral, which historically was used first in the study 
of these systems and which is still a very handy computational formula. 

Returning to (4), we have 

L,\.f(t)} = t, f(»T)c-" T \ (4) 

n— 

The inverse of the normal Laplace transform, with the simple change of 
variable t = nT and using p in place of the usual s, is 
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KnT) -i ( F(p)e ], " r dp, (10) 

where T, the contour of integration, encloses the poles of F{p), which is 
assumed rational. Substitution of (10) into (4) gives 

Uf{t)\ = £ e' nT ' £-. [ np)e pnT dp 

"=" m 2 *-> Jv (11) 

= o^Z [F{p)e l *- UT dp. 

Since the series (4) and the integral are both uniformly convergent, the 
order of summation and integration may be interchanged: 

Llf(t)) m±.[F(p)t, '■ U "" ) " r dp. (12) 

lirj Jr n =o 

If | e (p ~ t)T | < 1, we may sum the geometric series in (12) to 

Y „<"-">" r = * (13) 



ii =0 



1 _ e (l>s)T' 



where the condition that this summation be valid is that Re p < Re s. 
Thus, (12) becomes 

ujwi-r-./r-^U*- cm) 

The contour of integration, T, is the usual one parallel to the imaginary 
axis extending from — j« to -\-jac to include all possible poles. However, 
since we do not wish to exclude functions which do not vanish as s(or p) 
approaches infinity [in particular, F(p) = constant] the contour is closed 
by an infinite semicircle to the left from +j« to — joo . The requirement 
that Re /; < Re s is equivalent to stating that V shall include the poles 
fo F(p) but exclude the poles of 



1 _ e (/>-*> t • 

If the inverse Laplace transform of 

1 



F(s) = 



1 - c-» T 



be interpreted as the sum of a sequence of unit impulses, 8 T (l), a distance 
T apart, then the amplitude modulation of 8 T (t) by some function /(/) 
may be found, from the complex convolution formula of Laplace trans- 
form theory, to be: 
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L[f(t)8 r (0] = -L f — i 
27T/ Jr 1 — 



F(p) 



=Zi d V, 



c (p-a)T 



which is, of course, the Laplace-Stieltjes transform of f(t) as given by 
(14). Hence, the Laplace-Stieltjes transform is formally equivalent to 
the results of impulse modulation in the sense that both lead to the same 
transform. However, the definition of the Laplace-Stieltjes transform 
is rigorous and it directly relates discrete and continuous systems. In- 
tuitively, one would expect that, as the interval between samples ap- 
proaches zero, the Laplace-Stieltjes transform should approach the 
Laplace; i.e., the discrete system should look more like the continuous. 
This follows from the definition of the Laplace-Stieltjes transform. 
Note that, for the a(t) of Fig. 1, 

lim Ta(l) = t. 

Thus, 

lim TLlf(l)} = lim f f(t)r' sl dTa(t) 
r—o r-»o Jo 

= f f{t)c-* l tU = L[f(t)}, 

which is the desired relation, the Laplace-Stieltjes transformation ap- 
proaching the Laplace in the same manner as the staircase distribution 
function Ta(t) approaches the straight line /. 

Equation (14) is a very useful computational tool and can be used to 
prepare a table of Laplace-Stieltjes transforms from the common tables 
of Laplace transforms. Some elementary functions are given in Table I, 
where a direct comparison between the two transforms can be made, 
the Laplace-Stieltjes transform being written in the c* T form. Such a 
comparison is interesting, but the relationship between the two trans- 
forms can be better shown by a closer examination of the transforms of 
some elementary functions. Consider first /(/) = e~ al : 



1 



Lie- 1 ) = — '— = F(a), (15) 

S + a 

Ue- at ) = -A—-, = F(e' T ). (16) 



,,sT n—aT 



The single pole at .s = — o of /• , (.s) is shown in Pig. 2(a) in the usual man- 
ner. However, F(e sT ) has an infinite number of singularities occurring at 
s = —a ± 2irn/T(n = 0, 1, 2, • • • ). Thus, the effect of sampling is to 
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multiply the single real pole of the L transform into an infinite number of 
complex poles as shown in Fig. 2(b). The case of imaginary roots in the s 
plane is even more interesting. Let f(t) = sin fit 



L[sin 0t\ = 



Ljsin |8/] = 



s- + 13- 



5 = ns), 



e T sin &T 



„i*T _ Oo*T 



2e* T cos pT + 1' 



(17) 



(IS) 



The singularities of F(s) are shown in Fig. 3(a) and those of F(e"') in 
Fig. 3(b). Here, the sampling multiplies the original pair of poles into 
an infinite number of such pairs. The center of each pair is separated from 
the next by a distance of 2ir T. The distance or period T may be identi- 
fied with a radian sampling frequency w, = 2ir T. From Fig. !i(b) it is 
apparent that, if ^ w, 2, the pairs of poles overlap each other and form 
a new configuration which is indistinguishable from a configuration 
resulting from some function of radian frequency less than oj.,/2. This 
result also follows from Shannon's sampling theorem which, in effect, 
states that, if the original signal is to be recovered after sampling, then 
the sampling frequency must be greater than twice the highest frequency 
sampled. 

If the L transformation of a function has a singularity in the right half 
of the s plane, the L, transformation will have an infinite number in the 
right half plane. 



Table 1 



/(I) 


F<.) 


fdiT) 


FU'T) 


</(/) 


1 
8 


ii(nT) 


e .T 


,,-r _ ! 




1 


c -*»r 


e- T 




8 + <» 


e .T _ C -*T 





.s-s + 0* 


sin 0ii T 


e' T sin 0T 


sili 01 


e 2.T _ 2c" T cos pT + 1 


f <" sin 01 





p-anT s ; n pnT 


e i,+a)T s i n 1371 


(S + «)-' + 0-- 


e 2(.+<»»r _ 2<! <«"»>7' cos 07' + 1 


(■-■•' cos 01 


8 


,,-anT co g fj,,']' 


( e <«+.»r _ cos 0T)ci' +a,T 


(8 + a)- + 0- 


e i(ri a )r _ 2*<" + <" 7, cos 07' + 1 
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Fig. 2— (a) Single S plane pole of F(s) = l/(s 
complex S plane poles of F(e" T ) = e* T /(e? T — e~ aT ). 



a); (b) infinite number of 



Im s 



Im s 



J ft 
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Re s 



Fig. 3 — Pair of S plane poles of F(s) = 0/(s 2 + 2 ); (b) infinite number of S 
plane pole pairs of F(c" r ) = e" T sin 0T/(e 2 "T - 2e tT cos /3T + 1). 
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The change of variable, z = e sT , which introduces the z transform, 
simplifies the above mappings. It is well known that the transformation, 
z = e aT , maps the imaginary axis of the complex s plane into the unit 
circle about the origin in the z plane. The left half of the s plane maps 
wholely within the unit circle, the right half maps exterior to it. Since 
the z plane mapping is repetitive for multiples of 2w/T, the infinite num- 
ber of roots and root-pairs in the s plane map into single roots and root- 
pairs in the z plane. 

V. APPLICATION OK THE L 8 TRANSFORMATION TO THE SOLUTION OF 
FINITE DIFFERENCE EQUATIONS 

Many of the concepts of sampling can be applied to the solution of 
linear finite difference equations, with constant coefficients. These equa- 
tions are simply linear combinations of sequences of numbers shifted 
forward and backward in time by integral multiples of some fixed inter- 
val. In the case of a digital computer operating in a control loop, the 
sequences are actually generated by sampling some continuous function 
of time. If the assumption is made that all the sequences in a finite dif- 
ference equation result from such sampling, then the L, transformation 
offers a very useful method for the operational solution of such an equa- 
tion. The resulting solution is the "smooth" curve which, when sampled, 
will give the sequence satisfying the difference equation. 

A finite difference may be defined in either of two ways. One could be 
called a backward difference, defined as 

h[y(nT)\ = \!j(nT)\ - [y{nT - T)\. (19) 

That is, the backward difference is simply the difference of two sequences, 
one of which is the other shifted backward one interval in time. Since 
there is no possibility of ambiguity, the braces may be omitted and (19) 
written in the more usual form 

A b //„r = y„T — U„t-t ■ (19) 

In similar fashion, the forward difference may be written as: 

A/V/„T = UnT+T - U„T • (20) 

Higher differences of course are formed by taking "differences of differ- 
ences." That is, 

Ak"y,,T = ±b"~ l !J„T — Ab"~ 1 HnT-T 

or 

A/// H r = A/" -l //«r+r — A/" _1 ?/«r • 
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In order to eliminate possible confusion between differences, a difference 
equation can always be expanded and written in the ordinate form : 

n "t 

2 bjtjnT+jT + 23 CliUnT-iT = X nT . 

i=i i'=o 

If the assumption is made that the sequences [y(nT)}, \x{nT)) re- 
sult from a sampling of some continuous function which has an L s trans- 
form, then the finite difference equations of the above form are readily 
solved by the L s transformation. This follows from the results of Prop- 
erty I of the L„ transform, which is proved in the Appendix. The property 
is repeated below without proof. 

Property I: If L„[f(t)] — F{z) and a is a nonnegative integer, then: 

L a [f{t - aT)\ = L„[{f(nT - aT)}] 



= z- a \F{z) + Hf(-mT)z' 

and 

LJLftt + aT)] = LJLlftnT + aT)\] 

F(z) - £ f(?nT)z-' 



(21) 



= z 



(22) 



Application of (21) to the backward difference equation (14) leads to 

L 3 [A b Y nT ] = i^-i Y(z) - y(-T) (23) 

z 

and to the forward difference 

L,.[AAr] = (z - l)Y(z) - zy(0). (24) 

The terms y(0) in (23) and y{ — T) in (24) are the usual initial condi- 
tions, and allow the specification of arbitrary boundary conditions in a 
manner completely analogous to the insertion of initial conditions in the 
solution of differential equations. However, for simplicity, zero initial 
conditions will be assumed in the problem below. 

As an example, the above can be applied to the simultaneous difference 
equations : 

0.25ty„_i + w» = %n , 

(25) 
-1.0//,, - 0.25//„_i + w„ + 0.5«;„_! = 0, 

with T taken as unity for convenience. Letting 
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Y(z) = LXuX 

W{z) = L,[w n ], 

X(z) = L s [.r„], 

and performing the indicated operation on (25) leads to: 

Q.2oY(z) + zW(z) = zX(z), 

-{z + 0.25) Y(z) + (z + 0.5)W(z) = 0, 

whence 



(20) 



= z(z + 02o)X(z) (2?) 

K) z- + O.oz + 0.125 

As a "test function" let x(n) be the unit sample defined as unity for 
n = and zero elsewhere. The L„ transform of .i„ is then L,[x„] = 1, 
and it follows that 

WW = Z{Z + °- 25) • (28) 

WW [(2 + 0.25)2+ .0625]' 

By the application of (9), 

w„ = 4(0.35)" +1 sin (n + 1) ^ + (0.35)" sin n ^ . (29) 

In a similar manner w„(or y„) may be determined for any function x„ 
which is L t - transformable. 

In electrical network theory based upon the Laplace transform, the 
weighting function of a network is the time response of that network to a 
unit impulse. The analogy to the above response of a difference equation 
to a unit sample is clear. The fact that the difference equations could 
describe the operation of a digital computer performing linear operations 
in real time lends physical reality to the analogy and introduces the con- 
cept of a digital network. The designer of a servomechanism or feedback 
amplifier is concerned that his device shall be stable. The designer of a 
program for a digital computer is likewise concerned that his machine 
behave in a stable manner. Here, stability is defined in the sense of the 
electronic network designer that, for any bounded input, the output shall 
not continually increase. This is more elegantly and precisely stated in 
terms of the complex s plane. Here, the criterion for stability is that the 
characteristic equation have no roots in the right half of the 8 plane. 
Since the defining transformation z = e sT of the z plane maps the right 
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Fig. 4 — Nyquist plot of z 1 + 0.5* + 0.125. 



half of the s plane onto the exterior of the unit circle, the criterion for 
stability of finite difference equations is that the characteristic equation 
shall have no roots exterior to the unit circle in the z plane. Intuitively, 
it would seem that Xyquist's criterion could be applied in one form or 
another directly in the z plane to determine stability. However, the ap- 
plication is not as attractive as in the continuous case. The finite differ- 
ence equations (25) will be used to illustrate the point. As above, the 
roots of the denominator of (27), the L„ transformation of (25), determine 
the stability of the difference equations. Since all roots do indeed lie 
within the unit circle, as z takes on values along the unit circle in the 
positive soisc, the plot of z~ + 0.52 + 0.125 will encircle the origin in 
the positive sense a number of times equal to the number of zeros of the 
polynomial which lie within the unit circle (two in this case). Fig. 4 is 
such a plot and it is readily apparent that there are two encirclements. 
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In principle, the method can be extended to a polynomial of any order. 
As a practical matter, counting the number of encirclements from a poly- 
nomial of high order without making a mistake would be difficult. How- 
ever, this is not the most important point. One of the great attractions 
of the Nyquist criterion to the practical designer of feedback devices is 
that not only does it determine stability but it also indicates at a glance 
the margins against instability. To know the allowable variation in the 
gain of an amplifier, as a specific example, is of great value to the de- 
signer (and manufacturer) of a feedback amplifier. Unfortunately, in 
the discrete case no such information is apparent and the indication is 
simply one of "go" or "no go". 

If the change of variable w = z~ ] = e~ sT is made, the right half of the 
s plane is now mapped into the interior of the unit circle, and hence the 
criterion for stability becomes the more usual one of not enclosing the 
origin. Illustrating, the characteristic equation of (27) becomes; 

io~ 2 + 0.5uT' + 0.125. (30) 

As w takes on values on the unit circle in the positive sense, (30) has 
exactly the same values as shown in Fig. 4 with the exception that the 
curve now encloses the origin in the negative sense and, hence, stability 
depends also upon the sense of enclosure. The difficulties above will 
usually require that stability analysis still be made in the s plane. How- 
ever, some advantage can be taken of the angle preserving properties 
of the change of variable s = \/T In z. 

VI. CONCLUSION 

The Laplace-Stieltjes derivation of the z transform is straightforward 
and rigorous. As a mathematical model of the sampling process it has 
the advantage of also describing some of the operations possible within 
the computer. In particular, since it can be used as a basis for the opera- 
tional solution of linear finite difference equations it closely relates the 
solution of discrete linear systems and continuous linear systems. It 
is of considerable advantage to be able to apply the methods of network 
analysis to this type of computer operation. The Laplace-Stieltjes trans- 
formation forms the connection in a very clear manner. 
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APPENDIX 

PROPERTIES OF THE L s TRANSFORMATION 

The utility of the L s transform is increased by various properties per- 
taining to its use. The more important of these properties are derived 
below with some discussion of the area of application. Such discussion 
must of necessity be brief. 

Real Translation 

This property is the basis upon which the operational solution of linear 
finite difference equations is based. 
Property I : If L„[f(t)] = F(z) and a is a nonnegalive integer, then: 



F{z) + E/C-mTV"], 

m-1 J 

F(z) - Z f(mT)z-'"] . 



L.W - aT)] = z - 

and 

L s [f{t + aT)] = z a 

The proof of the first part follows: 
By definition 

U/Xr)] = f j{r)e- ST da{r) = F(e' T ). 
Jo 

Dividing the range of integration into two parts: 

F(e sr ) = r f(r)e " Mr) - f Mf" Mr). 

J-aT J-aT 

We now let t = t — aT in the first integral on the right-hand side : 

F( e ° T ) = f" f( ( - a T)e~ BU - aT) da(l - aT) - f /(r)e" ST da(r). (31) 

From Fig. 1 it is apparent that 

da(t - aT) = da(t). 
Hence, 

F(e sT ) = e aTs I f(t - aT)c7 st da(t) - [ Me'" Mr). (32) 



THE Z TRANSFORMATION 193 

But the second integral may be evaluated as: 

[ f( T )e- T clair) = J2f(-mr)e mTa . (33) 

J-aT m=l 

Substitution of (33) into (32) and rearrangement leads to 

<,-"• [~F(e ,r ) + Zf(-mT)e mT8 ~\ = £ f(t - aT)e~ at dait). 

But the right-hand side of the above is by definition L s [f(t - aT)]. 
Making the usual substitution of z = e T now leads to the desired result: 

L,[f(t - aT)] = z~ a If(z) + EjT-mTV'] . (34) 

We note in particular that, for a = 1, 

L.W -<&)] = z~\Hz) +K-T)}. 

For proof of the second part, the range of integration is divided into 
two parts: 

L.[/(t)] = f f(r)e" 8T da(r) + [ " /(r)r^(r). 

Letting r = t + aT in the first integral on the right-hand side, we have 

F ( e > T ) = f f( t + a7>"' < ' +a7 ' , dait + CIT) + [' fi^e'" da(r). 

Jo- J o 



Rearranging and substituting dait + aT) = da(t), 
e aT Tw r )- t" ' + fir)rr Bt dait) 



= f /(/ + «r)r r f/a(/). (35) 



The integral on the left is obviously 

f" F+ Me-" Mr) = Z f(»>T)c-' lT \ 

Jq m= 

and substitution into (35) gives 

e °r* [p^ _ g f( mT ) c - 

whence 

Wit + aT)] - z a \Fiz) - E /(m7V 

L PB-0 

For the special case of a = 1, we have 



= L a \fit + aT)], 



(3(i) 
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L.W + T)\ = z[F(z) - /(()) | = L a \f(t + 201- (37) 

Finite Differences 

An immediate consequence of the real translation property is that 
for finite differences. If finite differences are defined as in the text, we 
have: 

Property II: If the sequence {f(nT)\ resulting from the sampling of the 
continuous function f(t) has the L s transform F(z), then: 

L s [A b {f(nT)}} = L a \A b f nT ] = llli F(z) - f(~T), 

z 

L.(A// nT ) = (z - l)F(z) - 2/(0). 
By definition, 

L s (A 6 /„ r ) = L,(f» T - f„T-T). 

By linearity of the transform and (3-4), 

L^Ai/.r) = Fit) - z-\F{z) +/(0)] 

"( fL r^) F(f) _/( ~ n 

Again, by definition, 

L s (A//, ( r) = L s {f llT +T - f„r) 
and, by linearity of the transform and (37), 

LAArLr) = z\F(z) - ./'(()) | - F(z) - (z - l)F(z) - zf(0). (39) 

Finite Summation 

As integration is the inverse operation of differentiation, summation 
can be considered the inverse operation of taking differences. This is 
demonstrated more clearly in the property below. The process of finite 
summation is best demonstrated in the case of a computer operating 
in real time. The computer samples a function which is continuous and 
well defined at the sampling instants. At each sample the computer 
adds that sample to the sum of all preceding samples. If the result of 
this operation is a sequence \g(?it)}, we have as the value of gird) at any 
time nT: 

9(nT) = E/(fcT), 
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where {f(kT)}k is the sequence of sampled inputs. The analogy to inte- 
gration with respect to time is clear. The L 3 transform of such a summa- 
tion is given below. 

Property III: If [g(nT)\ is an infinite sequence such that each value of 
it is given by g(nT) = £"= f(kT), then G(z) = z/(z - l)F(z), where 
F(z) is the L t transform of the sequence \f{nT) } . 

By definition : 

G{z) = t, 9(nT)z-. 

n=0 

Substitution of the value for g(nT) gives 

„=o [_fc-o J 

Since we are dealing with uniformly convergent series, the order of 
summations may be interchanged, provided a suitable change in the 
limits is made in a manner equivalent to the change in limits when the 
order of integration is interchanged in double integration. Thus, 

g(z) = E/(at) Er, 

which may be written as 

g« = L/(at)«-* Er. 

However, the series in n may be summed as 1/(1 — z~ ), and hence 
G(z) = Zf(kT)z- k \-^- I, 

A=U \_Z — l_ 

and hence 

G(z) = -1- F(z). (40) 

2—1 

Complex Multiplication 

The superposition property of electrical networks is a very elegant 
and useful result of their linearity. For continuous linear networks, super- 
position is most concisely represented as a convolution integral, which 
has a particularly important Laplace transform. The same ideas also 
apply in the discrete case with the integral replaced by a summation. 

Property IV: Iff(t) and w(t) have the L, transforms F(z) and \Y(z), then: 

W(z)F(z) = L,\f:uikT)f(nT - IT) 
L*-° 

By definition, 



196 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1959 

W(z) = SwCMV*. 

A-=0 

W(z)F(z) = J2 w(kT)z" k F(z), 

A=0 

but 

z~ k F(z) = LJLf{t - kT)} - L a {[/(nT - kT)]}. 
Hence, 

TT(«) W = E w(kT)L a \f(t - kT)] = E L 8 [w(kT)f(t - kT)} 

A=0 fc=0 

(41) 

-iii[g«(*r)A< - *d]. 

At time £ = wT we have 

IF(i)F« = L. [g w(kT)f(nT - kT)~] 
but /(*) = 0; J < and therefore 

JF(z)F(z) = L, |"E w(kT)f(nT - fcT)l . (42) 

ScaZe Change 

Properly V:IfL.\f(t)] = F(z), then L.[e" a '/(«)] = F(kz), where k = e +aT . 
From definition, 

L.[e-°7(0]= C f(t)e- i ' Ta)T da(t) 

= E i^nDle-VT, 

n = 

Ue* T ffl] - Zf(nT)[kz}- n = F(kz). (43) 

71 = 

REFERENCES 

1. Ragazzini, J. R. and Zadeh, L. A., The Analysis of Sampled Data Systems, 

Trans. A.I.E.E., 71, Part II, November 1952", pp. 225-234. 

2. Barker, R. H., The Pulse Transfer Function and Its Application to Sampling 

Servo-Systems, Institution of Electrical Engineers (British) Monograph No. 
43, July 1952. 

3. Truxall, J. G., Automatic Feedback Control Si/stem Synthesis, McGraw-Hill, 

New York, 1955. 

4. Jury, E. I., Analysis and Synthesis of Sampled Data Control Systems, Trans. 

A.I.E.E., 73, Part I, September 1954, pp. 332-346. 

5. Widder, D. V., Advanced Calculus, Prentice-Hall, New York, 1942. 

6. Gardner, M. F. and Barnes, J. L., Transients in Linear Systems, John Wiley & 

Sons, New York, 1942. 

7. Widder, D. V., The Laplace Transform, Princeton University Press, Princeton, 

N. J., 1946. 



